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Обсуждаются вопросы планирования эксперимента в интеллектуальных системах исследовательского 
проектирования морских судов. Особое внимание обращается на решение проблемы планирования 
эксперимента при реализации принципа конкуренции в условиях неопределенности и неполноты 
исходной информации. 

Ключевые слова: интеллектуальная система, планирование эксперимента, принцип конкуренции 


Оцезйоп$ оЁ р1апиих ехреттеп шт шеШоепсе зузетаз оЁтезеагсВ Чез1еие оЁР таппе $1р$ аге 41зсиззе4. Тве 
зресла! айепйоп а4@геззез оп Ше 4еслз1оп оЁ а ргоет оЁ р!апииз ехреттепи аё геапхайоп оЁ а сотрейбоп 
рипстр1е ш Ше сопаюоп$ оЁ ипсецапиу апд пшсотр!еепез$$ оЁ фе шша[ шРогтайоп. 

Кеу могд5$: пиеПесваа| зузет, 4ез1еп оЁРехреттеп($, пе рипстре оЁ сотреНвоп. 


Обговорюються питання планування експерименту в 1нтелектуальних системах дослдного проектування 
морських суден. Особлива увага звертаеться на розв’язання проблеми планування експерименту при 
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1 Концепция моделирования экстремальных ситуаций 

В задачах исследовательского проектирования 


При реализации интеллектуальных технологий в интеллектуальных системах (ИС) 
исследовательского проектирования на базе высокопроизводительных средств вычис- 
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лений учитываются общие принципы и структура информационной модели, обеспе- 
чивающей анализ и прогноз исследуемых ситуаций. Такая модель разрабатывается в 
рамках нечеткого логического базиса и формализованных методов оценки динамики 
взаимодействия судна в различных условиях эксплуатации, в том числе и в экстре- 
мальных ситуациях [1-15]. 

В настоящей статье рассматривается приложение разработанной концепции интер- 
претации текущих ситуаций в сложных динамических средах [9]. Объектом исследо- 
вания выбрана область практических приложений, связанная с контролем динамики 
судна на волнении на этапе исследовательского проектирования (рис. 1). 

Задача выбора альтернатив и принятия решений в нечеткой среде связана с по- 
строением варианта (выбора альтернатив), удовлетворяющего заданным ограничениям 
и являющегося способом достижения цели. В процессе решения задачи, представ- 
ленной кортежем [7], [11]: 


<А, \, 0 >, (1) 


где А — множество альтернатив; \ — характеристики внешней среды; $ — система 
предпочтений; Ц — действие, необходимо выполнить определенную операцию над мно- 
жеством альтернатив, установив наиболее предпочтительную из них. 

Решение задачи (1) ведется с использованием нечетких графов и нечетких отно- 
шений [2]. Концептуальный базис нечеткой среды моделирования и планирования 
эксперимента основан на использовании парадигмы обработки информации в мульти- 
процессорной вычислительной среде [8] и достижений в области интеллектуальных 
технологий ХХ] века [7], [11]. Принципы преобразования информации в трудноформа- 
лизуемых средах сформулированы в работах [2], [9]. 


Преобразование информации 
моделирование -— прогнозирование -— принятие решений 


Интеграция Создание Формирование 
идей и НОВОЙ задач 


методов техники и стратегической 
исследования технологий инициативы 


Рисунок 1 — Информационная модель создания новой техники и технологий 
в ИС исследовательского проектирования 


Формализованное ядро системы интеллектуальной поддержки процессов построе- 
ния и использования моделей знаний при анализе текущих ситуаций реализуется в 
рамках нечеткого логического базиса. Фундаментальной основой такой интерпретации 
является концепция нечетких целей и ограничений [2], [7], [11]. 

Повышение эффективности функционирования процедурной компоненты базы 
знаний интеллектуальной системы (ИС) достигается за счет использования принципа 
конкуренции и формализации процесса обработки нечеткой информации [1] в высоко- 
производительных вычислительных средах. Другими принципами повышения эффек- 
тивности ИС новых поколений являются принцип открытости, принцип сложности и 
принцип нелинейной самоорганизации [9]. Реализация указанных принципов осуще- 
ствляется в рамках концепции мягких вычислений [15], интегрирующей нечеткую 
логику [1], ИНС [7] и генетический алгоритм [4]. 
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2 Реализации принципа конкуренции 
при анализе альтернативных решений 
в ИС исследовательского проектирования 


Интерпретация текущих ситуаций в сложной динамической среде осуществляется 
в рамках принципа конкуренции [8], обеспечивающего выбор предпочтительной вы- 
числительной технологии обработки информации методами теории катастроф [9] на 
основе концепции мягких вычислений [15]. 

Рассмотрим прикладную теорию: концепцию принципа конкуренции (СотренНоп 
Рипсре СопсерЕ- СРС) как систему вида [8], [9]: 


СРС > 5: $ = <В,Р, >, (2) 


где $ — система, отображающая теорию СРС; К -— классы всех операций и матема- 
тических зависимостей, описываемых СРС; Р — совокупность операций на множествах К; 
О — множества отношений между элементами класса К. 

На основе отображения (2) представим компьютерную интерпретацию теории 
СРС как 


КСРС)у5*: 5* = <В*.Р*,.О*>, (3) 


Здесь символ * указывает на различия между указанными множествами, опреде- 
ляемые особенностями компьютерной интерпретации. Условия сохранения целостности 
информации о динамике взаимодействия обеспечиваются путем реализации парадигмы 
обработки информации в мультипроцессорной вычислительной среде [8]. Необходимая 
производительность достигается с помощью процессоров, которые являются узлами 
вычислительной сети единой распределенной системы. При реализации принципа 
конкуренции используется бортовой многопроцессорный комплекс стандартной архи- 
тектуры [2], [7], ПП. 

Выбор предпочтительного решения в рамках принципа конкуренции осуществ- 
ляется подсистемой поддержки принятия решений. Постановка задачи (1) реализуется 
на основе подхода Беллмана -— Заде и метода иерархий (парных сравнений) Саати [2]. 

Подход Беллмана - Заде. Рассмотрим отображение ф:Х—>У, значениями которого 
являются элементы множеств Х и У. Это отображение понимается как реакция системы 
на входные воздействия хеХ или как некоторые оценки выбора соответствующих 
альтернатив. Нечеткая цель задается в виде нечеткого подмножества универсального 
множества реакций У, т.е. в виде функции принадлежности (ФП) ис:У>[0,1]. Решение 
задачи состоит в достижении цели при заданных нечетких ограничениях. Пусть неко- 
торая альтернатива х обеспечивает достижение цели со степенью ис(х) и удовлетворяет 
ограничениям со степенью ис(х). Тогда степень принадлежности альтернативы х реше- 
нию задачи равна минимальному из этих чисел, т.е. ФП решения: 


но(х) = п {ис(х), исх}. (4) 


При наличии нескольких целей И отраничений, а также с учетом их важности А 
и у нечеткое решение описывается соответствующими ФП, а неопределенность выбора 
альтернативы с максимальной степенью принадлежности нечеткому решению 19) опре- 
деляется условием: 


таххх ир(х) = тах, к шщ{ис(х), ис). (5) 
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Метод анализа иерархий. Реализация этого метода связана с выделением прио- 
ритетов (весов) признаков в целях выбора наилучших из них с использованием алгеб- 
раической теории матриц и экспертных процедур: 


Ут =Л дл; 

Ли °° ЛПы 81 (6) 
пл ЁРЫ:| 

Ли ых Л т бт 


где \/ — обратно симметричная матрица значений парных сравнений признаков 
относительно данного атрибута; л — нормированный вектор весов признаков; А, шах — 
наибольшее собственное значение матрицы \; П- результат определения глобальных 
приоритетов признаков Ш,..., Им; М — число признаков; лу (ЕТ,...п, ]=1,...101) — отно- 
сительный вес 1-го признака по ]-му атрибуту; 2; — относительный вес ]-го атрибута. 

На множестве альтернатив каждому критерию ставятся в соответствие лингви- 
стические переменные (ЛП) с терм-множеством Т;={Ти,...,Ги}, которые характеризуют 
предпочтения конструктора. Базовым множеством ЛП является допустимое множество 
значений критерия. Оценка эффективности альтернативы по отдельному критерию 
осуществляется с помощью дополнительной ЛП с терм-множеством $={51,...,3, мощ- 
ность которого совпадает с Т;. Введенная таким образом ЛП инвариантна по отно- 
шению к каждому из критериев СК, а ее базовое множество принимает значения на 
отрезке [0,1]. Между двумя ЛИ существует причинно-следственная связь в виде нечет- 
кого импликативного отношения: 


ивиСьУ): Х,хУ-> [0,1], х:Хьуеу, (7) 


где Х; и У - базовые множества значений 1-го критерия и оценки альтернатив. 

При функционировании ИС вычислительный эксперимент выполняется в виде 
следующих этапов: формирование альтернатив (сценариев экстремальных ситуаций); 
оптимизация и тестирование вариантов; анализ альтернатив и выбор предпочтитель- 
ного решения. 

На основании подхода Беллмана - Заде и метода анализа иерархий разработан 
алгоритм и программное средство моделирования и интерпретации экстремальных 
ситуаций при реализации парадигмы обработки информации в мультипроцессорной 
вычислительной среде в ИС исследовательского проектирования. Наряду с подходом 
Беллмана — Заде и методом анализа иерархий Т. Саати при интерпретации решений в 
системе интеллектуальной поддержки конструктора могут использоваться и другие моде- 
ли. Среди них следует выделить функцию выбора, модель выбора и метод эталонов [11]. 


3 Особенности планирования эксперимента ИС 
исследовательского проектирования 


При разработке сценариев взаимодействия судна с внешней средой широко исполь- 
зуются методы планирования эксперимента. Содержательный анализ проблем модели- 
рования, определяемых соотношением между реальным экспериментом, модельным 
экспериментом и теорией приводит к поиску стратегий моделирования с учетом осо- 
бенностей исходной информации и производительности вычислительных средств. Фор- 
мальные процедуры, лежащие в основе функционирования ИС исследовательского 
проектирования, предполагают широкое использование измерительной информации 
и методов математического моделирования. Решение прикладных задач анализа и 
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прогноза ведется с привлечением теории планирования эксперимента. Информация, 
получаемая в процессе нормальной эксплуатации ДО, представляет собой данные пас- 
сивного эксперимента. Активный эксперимент используется только на этапе тестиро- 
вания и адаптации ИС. При использовании данных активного эксперимента проблема 
выбора наилучшего математического описания может быть решена классическими 
методами планирования. В случае пассивного эксперимента эта проблема требует 
специальных подходов, учитывающих адекватность описания, дисперсию оценок и 
вычислительные трудности из-за вырождения информационной матрицы [5]. Анализ 
подходов к планированию эксперимента и поиску эффективных процедур обработки 
информации в сложных ИС реального времени содержится в работах [2], [7], [9], [11]. 

Возникающие при практическом использовании ИС исследовательского проек- 
тирования задачи планирования эксперимента часто выходят за рамки классического 
подхода и требуют специального обобщения с учетом особенностей рассматривае- 
мых проблем. 

Рассмотрим определения, связанные с построением планов эксперимента при 
реализации моделей интерпретации динамических ситуаций [3], [5], [9]. 

Определение 1. Под планом эксперимента & будем понимать вероятностную меру 
на области планирования И с конечным носителем зир & = {йа/, &(й)==0\. Обозначим 
через Н! конечномерное подпространство некоторого пространства Н (Н‚<Н), которому 
принадлежит функция регрессии, и будем считать, что по условиям эксперимента могут 
быть использованы только планы из допустимого множества 


Л := {& Ям: сага (зар &) < п}, (8) 


где и — число, характеризующее ограниченность ресурсов. 

Определение 2. Выбор пространства оценивания Н!=НЕ и оператора оценивания 
5=5; обеспечивает при фиксированном плане 6х нахождение наилучшего в метрике 
пространства Н приближения Е\ с Н: к произвольному элементу \ из Н. При этом 


пространство оценивания Не должно обеспечивать построение оценки т неизвестного 


элемента п «Н. 

Определение 3. Задача выбора процедуры л=(&, Ну, 5) восстановления 1 из Н 
на множестве допустимых процедур п является двухкритериальной. В этих условиях 
поиск оптимальной процедуры л удобно вести на основе решения задачи оптимизации 
с приоритетом, учитывающим систематическую ошибку. 


к = Аго шЕ В. (т), (9) 
где у — вероятностная мера для осреднения систематической ошибки В(&, Н1, 5) 


В,(&Н,5):= | В(&, Н, $) (т); 


В(ЕН,5):= ЧР (п, Ей): = |п- Ей 


2 
Н ; 


5 — оператор, с помощью которого производится оценка п =1(у) неизвестного 
элемента | ‹Ы. 

Определение 4. Для характеристики случайной ошибки используется функционал 
Ф(Е,Н\,5) от корреляционного оператора оценки т (определитель де Р(\), след ря) 
или другие критерии), и тогда 

п = АтешЕФ(п), (10) 


ЕП 
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где П =П- множество процедур п’, являющихся решением задачи (9). 
Определение 5. Задача (8) на множестве процедур 


=* ен: {Е в /»:Ё = АВШЕВ, (т. (11) 


рассматривается в предположении, что носитель плана 5ир & определен однозначно 
из решения задачи (10), но имеется свобода в выборе весов наблюдений Р; (} = 1, т). 
Это позволяет минимизировать случайную ошибку при следующих условиях 


* 


* 
я (Р):=(@® у... [С ыы 7 =Аго ШЁ Фен, (12) 


лед Р®'; Р; > 0, ]=1 т, о }; РЕчь) — корреляционный оператор, опреде- 


ляемый как ЛД = М7 (МЕ- информационный оператор). 

В прикладных задачах планирования эксперимента важное значение приобретает 
интерпретация таких свойств, как ортогональность и ротатабельность. При равноточ- 
ных измерениях эти свойства приводят к независимости коэффициентов регрессии 
(их одинаковой точности), что имеет принципиальное значение в поисковых процеду- 
рах Бокса-Уилсона при экстремальном планировании. 

Выбор оптимальных условий эксперимента. Задача выбора оптимальных усло- 
вий эксперимента в ИС определяет надежную оценку характеристик ДО и параметров 
внешней среды и связана с построением нормированного дискретного плана для ДО, 
развивающегося во времени и пространстве. 

Определение 6. Стратегия выбора оптимальных условий эксперимента харак- 
теризуется тем, что объект функционирует в режиме нормальной эксплуатации. ИС на 
основе анализа ситуации выбирает моменты времени и координаты точек, в которых 
следует производить измерения. Выполнив серию опытов при некоторых фиксиро- 
ванных значениях исследуемого фактора в различные моменты времени и имея модель 
системы, можно подсчитать нормированную информационную матрицу 


М(=)=УР,М(6,) (13) 
У=1 
для дискретного плана эксперимента 


а. - 
А = (14) 
Е _ | Уве Р. 50) 


п 


где Р‚ — веса наблюдений. 
Определение 7. Синтез Р-оптимальных планов измерений осуществляется на 
основе итерационной процедуры с использованием функции 


5((,)=зар[М" (&М(0]. (15) 


Определение 8. Особым случаем идентификации динамических систем является 
планирование измерений в частотной области. Управление процессом измерений на 
основании методов планирования эксперимента позволяет свести задачу к поиску 
дискретного оптимального плана 


а 1»... (16) 
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с нормированной информационной матрицей 


+. т. 
М(О = ро (ФТ в ‚ЭР, + 1) 


п 
1 Уфо ФТ юр, 
2-1 
где в, — частота спектрального разложения стационарной случайной функции 
(входного процесса); Р, — ординаты непрерывного спектра этой функции; звездочкой (*) 
помечены комплексно сопряженные частотные характеристики. 
Интервал, на котором определена вероятностная мера, порожденная нормиро- 
ванной спектральной плотностью, характеризуется выражением 


0 = [-бь, +04. (18) 


где [-ск, +ож] — диапазон частот, определяющий полосу пропускания. 

Спектр плана в этом случае представляет собой совокупность значений частот, 
а дисперсии гармоник являются весами. 

Практическое применение методов планирования измерительного эксперимента 
реализовано в виде встроенной процедуры в системе автоматизации эксперимента 
на стадии исследовательского проектирования. Процедура реализует выбор оптималь- 
ных условий измерений, обеспечивающих надежную оценку характеристик ДО и пара- 
метров внешней среды, а также формализацию процедурной компоненты базы знаний 
ИС и соответствующего алгоритмического и программного обеспечения. Следующие 
этапы эффективной реализации измерительного процесса связаны с разработкой ин- 
формационной технологии оперативного контроля состояния ДО и прогноза его пове- 
дения в различных условиях эксплуатации, в том числе и в экстремальных ситуациях. 
Конкретное наполнение знаниями прикладной области экспертизы определяется осо- 
бенностями измерительного процесса и решаемых задач с учетом целей и наложенных 
ограничений. 


4 Интеллектуальные технологии альтернативных 
решений на основе методов планирования эксперимента 


Интеллектуальная поддержка планирования эксперимента в зависимости от осо- 
бенностей исследуемой ситуации позволяет давать полезные советы по принятию 
решений на неформализованных этапах процесса экспериментального исследования, 
включая идентификацию ДО, выбор модели взаимодействия, выбор критерия оптималь- 
ности плана эксперимента, процедуры рандомизации и схемы эксперимента. В качестве 
модели представления знаний используется семантическая сеть. Ядром системы интел- 
лектуальной поддержки является набор функциональных блоков, работающих со зна- 
ниями, хранящимися в семантической сети. В вершинах сети располагаются модели 
представления четких и нечетких знаний и необходимая информация для выполнения 
операций по планированию эксперимента. Компьютерная программа обеспечивает 
работу системы в виде последовательной активации вершин сети, содержащих списки 
ссылок на другие вершины. 

База знаний интеллектуальной поддержки планирования эксперимента содержит 
различные типы знаний, среди которых выделяют знания описательного характера 
(факторы, источники неоднородностей и планы экспериментов) и логические правила, 
обеспечивающие принятие решений на основе анализа посылок, а также знания, содер- 
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жащие результаты компьютерного моделирования различных ситуаций при принятии 
решений о выборе плана эксперимента с учетом экономических, статистических и инфор- 
мационных требований. Механизм логического вывода поддерживается стеком актив- 
ных вершин со списком номеров вершин, ожидающих очереди активизирования. База 
данных системы содержит данные, используемые в процессе логического вывода, и 
вспомогательные базы, включающие в себя каталог планов эксперимента, библиогра- 
фическую информационную систему и каталог теорем. Программа для генерации планов 
эксперимента хранится на диске в отдельных файлах. 

Моделирование текущих ситуаций обеспечивается путем решения задач взаимо- 
действия судна с внешней средой, представленных в виде различных математических 
моделей, в том числе и блочно-модульной структурой. При этом исследуемая задача 
в заданной ситуации описываются в виде спецификации модели рассматриваемой 
области знаний. По спецификации аппарат предметного уровня формирует прикладную 
модель и перевод ее в контексте формального уровня модели предметной области, с 
помощью которого формируется модель ситуации и алгоритмическая ее реализация 
в виде набора прикладных программ. 


5 Нечеткий логический базис при анализе 
конкурирующих вычислительных технологий 
методами планирования эксперимента 


Рассмотренные в настоящем разделе методы планирования эксперимента преду- 
сматривают возможности использования достижений этой теории в различных задачах 
поддержки принятия решений. Однако в экспериментах, связанных с интерпретацией 
текущих ситуаций, особенно при сложном взаимодействии судна с внешней средой, 
возникают определенные трудности. Именно поэтому новый подход к построению 
алгоритмов принятия решений должен предусматривать построение плана экспери- 
мента в рамках принципа конкуренции, позволяющего выбирать предпочтительную 
вычислительную технологию. Поскольку принцип конкуренции основан на развитии 
концепции мягких вычислений, при адаптации стандартных планов эксперимента при- 
ходится учитывать особенности нечетких и нейросетевых моделей, а также алгоритмы 
эволюционного моделирования. 

Приведем примеры типичных логических правил, используемых при реализации 
принципа конкуренции. 

1. «ЕСЛИ число исследуемых факторов не превышает п=5 И все факторы 
представляют собой однородные экспериментальные единицы И данные эксперимента 
представлены нелинейной моделью, ТО рекомендуется использовать планы Коно или 
центральные композиционные планы». 

2. «ЕСЛИ опыты проводятся на неоднородных экспериментальных единицах 
И источники неоднородностей можно группировать в подмножества-блоки, ТО реко- 
мендуются блочные рандомизированные планы». 

Механизм логического вывода поддерживается стеком активных вершин со спис- 
ком номеров вершин, ожидающих очереди активизирования. 

Исходная математическая модель. Выбор исходной математической модели, 
описывающей динамику взаимодействия в нечеткой среде при планировании экспери- 
мента, состоит в определении вида функции [1], [9]: 


У=ЁЕ(®), (1=1,...,хк). (19) 
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Геометрический образ выходной характеристики у рассматривают как функцию 
или поверхность отклика, а функцию ЁКх;) принимают в виде отрезков степенных 
рядов (алгебраических полиномов) различных степеней. Функции Ё(х) при ]=1,2....,К 
образуют матрицу известных коэффициентов Е, которая полностью описывается выбо- 
ром модели Ё(х). По данным п измерений величины у; для независимых переменных 
х; (1=1,2.....п) в принятой области планирования можно составить матрицу Е и вектор 
У результатов наблюдений: 


оо аа | У 
О В (20) 
О У, 


Матрица Е и вектор У полностью описывают информацию для определения 
неизвестных коэффициентов математической модели, интерпретируемой в рамках теории 
нечетких множеств. Выбор числа и условий эксперимента, обеспечивающих получение 
наилучшего результата, составляет цель планирования эксперимента в нечеткой среде. 
В качестве критерия оптимальности в задачах планирования нечеткой системы исполь- 
зуется О-оптимальность. Учет случайной и систематической погрешности произво- 
дится на основе обобщенного критерия оптимальности плана Х*, минимизирующего 
на множестве планов выражение 


р 2 


[Е У.7096, - (а, | 4, (21) 


[4 
в котором в качестве модели у 7. (х)0, может быть выбран полином меньшей 
1=1 
степени, чем в истинной модели. 

Проблемы теории планирования эксперимента при интерпретации текущих си- 
туаций в задачах принятия решений связаны со сложностью процессов взаимодей- 
ствия судна с внешней средой в условиях неопределенности и неполноты исходной 
информации. Процессы функционирования в условиях такого взаимодействия носят 
динамический характер, а измерения параметров определяются повышенными требо- 
ваниями к их точности. 

Оптимальная структура модели. Рассмотрим в качестве исходной математи- 
ческой модели зависимость [3], [5], [13]: 


У=хВ+еЕ; (22) 
2 

Е(&) = 0, соу(&) = о Г (23) 
где Х — (№хК) — матрица наблюдений регрессоров х;; В — вектор-столбец значений 
неизвестных коэффициентов размерности (^+1); = — вектор-столбец значений регрес- 
сионной ошибки размерности (^+1); У — вектор-столбец выходной переменной; Е — опе- 
ратор математического ожидания; /— единичная матрица; № — объем выборки; К — число 

регрессоров:; о? — дисперсия воспроизводимости. 


Точность прогноза модели (22) определяется дисперсией ошибки предсказания. 
Для (№+1) наблюдения по адекватной модели эта дисперсия равна 


—1 
АИ. 7 о? х,.д(Х7Х) ен 2 (24) 
где Хм. 1х — вектор-строка размерности (1х). 
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Для сопоставления дисперсий моделей с «недобором» и «перебором» регрес- 
сионных переменных относительно истинной зависимости выражение (24) преобра- 
зуется к виду: 


ВАУ, ] = ое С 


2 
т В, / 2 

Хуа я П + Ил? Аим 

/ к 


+1 


(25) 


где Ак — параметр нецентральности {-распределения с №-К степенями свободы 


= [2/5 [(н-00-# |; (26) 


с, - выборочная дисперсия переменной хк; А, - выборочный коэффициент мно- 


жественной корреляции переменной хк с остальными независимыми переменными 
хр, ..., Лк, ВЫЧИСЛенный по данным матрицы наблюдений Х. 
Дисперсия ошибки прогноза с «недобором» и «перебором» будет равна 


ВАУ, | — ее о _ 
Хунин + 1+ В, о т а) 


ВАУ), | = р а Х 


(27) 


Г 2/2 > ро 
Хун + П+В,/^, а = о 


Здесь Вк — значение коэффициента при А-й неучтенной переменной хх; хмик — зна- 
чение (№+1)-го наблюдения (-й переменной; 


Хи = Хмщьт Ак - регрессия переменной х, на остальные переменные 

К 
= АХ ХХХ — вектор-столбец коэффициентов размер- 
Хх = Я Ха? АЯ АА КА 

Ч 


ности (К-1)х1. 
На практике приходится избегать переусложнения модели за счет включения 
дополнительных регрессоров, так как 


РАУ,. |< ВАУ. (29) 


что следует из сопоставления (24) и (25). При этом равенство дисперсий дости- 
гается только при 


Хмнкы = Амьен. (30) 


Выбор наилучшего уравнения поверхности отклика может быть произведен также 
с использованием базисных функций [3], [5] и модели вида 


К 
и: У Вх № (и), (30 
К=1 


где Вк — К-й коэффициент регрессии; (и) — К-я базисная функция (В (и)=1); и = 
т 
(и, ..., им) — вектор факторов пассивного эксперимента; М — число факторов; К — общее 
число слагаемых. При этом факторы кодированы, а модель наблюдения имеет вид 


У=п +=. (32) 
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Системный анализ задачи выбора наилучшей полиномиальной регрессии позволяет 
сформулировать подход к определению степени тор, алгебраического многочлена Р„(с,х), 
достаточно надежно аппроксимирующего заданную экспериментальную зависимость 
[Ух]. Суть его состоит в дополнении классического метода (выбор в качестве тор, 
величины т, доставляющей минимальное значение сумме квадратов отклонений) 
описанием способа измерения экспериментального массива {у„} 


У = 8(С (Хи) + 2и), (33) 


где С(х„) - некоторая функция; &„ — случайные ошибки (п = 1,№). 

Особенности вычислительной технологии. Вычислительные трудности при реали- 
зации алгоритмов обработки информации в трудноформализуемых средах, особенно 
данных пассивного эксперимента, связаны с процедурой обращения матрицы Х"Х 
системы нормальных уравнений 


(ХХ)! В = ХТУ. (34) 
Как показывают результаты вычислений, обусловленность информационной мат- 
Т. 
рицы ХХ существенно хуже, чем матрицы Х. Если Лах И Ап — наибольшее и наимень- 
и 2 Т 
шее сингулярные числа матрицы Х, а Аки ЛА — матрицы ХХ, то при небольшом 
2 > й 
Ам число Ам может иметь Порядок, сравнимый с порядком погрешностеи. В результате 


матрица ХХ становится почти вырожденной, тогда как матрица Х таким свойством 
не обладает. При решении плохообусловленных задач используют методы псевдообра- 
щения, обеспечивающие работу вычислительных процедур в условиях вырожденности. 
Псевдообратная матрица Х’ существует для любой матрицы Х, и для любого вектора 


^ ь й о 
УВ = Х У с минимальной нормои среди всех векторов, минимизирующих сумму 
2 


квадратов отклонений |у — ХВ 


Сравнительный анализ экспериментальных данных с помощью различных под- 
ходов обращения матриц позволяет выделить прямой метод Гревилла, позволяющий 
работать непосредственно с матрицей Х. Эта матрица лучше обусловлена, чем матрица 
ХТХ. В результате существенно повышается устойчивость решения. При использовании 
нелинейных по параметрам моделей возникает проблема оценки коэффициентов рег- 
рессии. При некомпактности априорного множества параметров такая оценка может 
вообще не существовать (не достигается минимум суммы квадратов отклонений). Алго- 
ритмы минимизации суммы квадратов основаны на методе Ньютона -— Гаусса и сво- 
дятся к линейной аппроксимации функции регрессии. 

Планирование эксперимента в ИС исследовательского проектирования связано с 
выбором оптимальных условий, обеспечивающих организацию процедурной компонен- 
ты базы знаний и соответствующего алгоритмического и программного обеспечения. 
Конкретное наполнение знаниями прикладной области определяется структурой ис- 
следуемого процесса взаимодействия и номенклатурой решаемых задач с учетом 
поставленных целей и наложенных ограничений. 

Практические рекомендации. Приведем некоторые рекомендации по использо- 
ванию методов планирования эксперимента в ИС исследовательского проектирования. 
Численный пример оценки критериальных соотношений для различных способов фор- 
мирования регрессионной структуры приведен в табл. 1. Здесь О[ГРУх.1] — дисперсия 
ошибки; \/ — частота; К; — выборочный коэффициент множественной корреляции; $30 — 
остаточная сумма квадратов; | — полная адекватная модель; 2 — модель с «недобором»; 
3 — модель с «перебором». 
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Особенности планирования эксперимента при моделировании... 


Таблица 1 — Характеристики различных способов формирования 
регрессионной структуры 


Вид Рассматриваемые критерии 
модели 


Оценка | РУ в | % | 


Ги [027 [055 | 050 [046-107 | 


В зависимости от характера исследования эксперимент можно выполнить либо 
в виде последовательной серии опытов, либо случайным образом (рандомизированный 
эксперимент). Типичными планами эксперимента в задачах исследовательского проекти- 
рования являются матрицы плана Коно и центрального композиционного плана [6], [13], 
наиболее часто используемых в задачах интерпретации ситуаций в сложной динамиче- 
ской среде. Для построения оптимального плана эксперимента необходимо задаться 
видом модели, определить область варьирования переменных и выбрать число испы- 
таний п. Если полученное значение п кратно двум или четырем, то для построения 
оптимального плана можно использовать матрицы Адамара. 

Если помимо линейных эффектов исследуемая модель содержит взаимодействие 
факторов, то при построении плана эксперимента с целью сокращения числа измерений 
удобно использовать дробные факторные планы, а в отдельных ситуациях — неполно- 
блочные сбалансированные планы (квадраты Юдена). Если исследуемая функция неиз- 
вестна, то можно ограничиться применением стандартного плана, при котором последо- 
вательно варьируется каждая из независимых переменных при условии неизменности 
всех остальных. 


Заключение 


Для развития науки исследовательского проектирования характерно противо- 
речивое единство двух тенденций: дифференциации и интеграции. Дифференциация 
обусловлена различием в объектах и методах теоретических и экспериментальных ис- 
следований, спецификой используемых научных дисциплин и приложений. Тенденции 
к интеграции знаний проявляются в росте значения науки исследовательского проек- 
тирования как единой системы знаний, в распространении общих методов исследования 
и их математизации. Развитие интеллектуальных технологий в условиях неопределенно- 
сти и неполноты исходной информации требует эффективного использования методов 
планирования эксперимента в ИС исследовательского проектирования на базе высоко- 
производительных средств обработки информации в сложной динамической среде. 

Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры 
инновационной России» на 2009 — 2013 годы (ГК № П295 от 30.04.2010, ГК № П976 
от 27 мая 2010, ГК № 02.740.11.0837 от 11.06.2010). Результаты работ использованы 
в рамках проекта реализации Правительства РФ № 218 «О мерах государственной под- 
держки развития кооперации российских высших учебных заведений и организаций, 
реализующих комплексные проекты по созданию высокотехнологичного производства». 
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